NONLINEAR SCHRODINGER EQUATION WITH A WHITE-NOISE 
POTENTIAL: PHASE-SPACE APPROACH TO SPREAD AND SINGULARITY 



ALBERT C. FANNJIANG 

Abstract. We propose a phase-space formulation for the nonlinear Schrodinger equation with a 
white-noise potential in order to shed light on two issues: the rate of spread and the singularity 
formation in the average sense. Our main tools are the energy law and the variance identity. The 
method is completely elementary. 

For the problem of wave spread, we show that the ensemble-averaged dispersion in the critical 
or defocusing case follows the cubic-in-time law while in the supercritical and subcritical focusing 
cases the cubic law becomes an upper and lower bounds respectively. 

We have also found that in the critical and supercritical focusing cases the presence of a white- 
noise random potential results in different conditions for singularity-with-positive-probability from 
the homogeneous case but does not prevent singularity formation. We show that in the supercritical 
focusing case the ensemble-averaged self-interaction energy and the momentum variance can exceed 
any fixed level in a finite time with positive probability. 



1. Introduction 

High-intensity laser beams propagating in dielectric media such as optical fibers, films or air 
are important problems both from fundamental and practical perspectives. The physical process 
is nonlinear and the amplitude modulation is customarily described by a nonlinear Schrodinger 
equation with an additional inhomogeneous potential V representing the random impurities in the 
medium. 

Let z and x E M, d be the longitudinal and transverse coordinates of the wave beam. In the 
physical setting, the transverse dimension d may be 0, 1,2. The simplest non-dimensional form of 
the nonlinear Schrodinger (NLS) equation with a random potential V then reads 

(1) i— ^ + -A x ^ + 7- 1 5|*| 2fT * + 7" 1 V r o^ = 0, a>0. 

oz 2 

Similar equations also arise in many other contexts such as in plasma physics and the Bose-Einstein 
condensation, see ^3] and references therein. Here we have written the equation in the comoving 
reference frame at the group velocity and non-dimensionalized the equation with the longitudinal 
and transverse reference lengths L z , L x ; 7 = L z j (kL%.) is the Fresnel number; g is nonlinear coupling 
coefficient with g > representing the self-focusing case and g < the self-defocusing case. For a 
nonlinear Kerr medium a = 1 leading to the cubic NLS equation. We shall let a be an arbitrary 
positive constant. We shall consider white- noise-in-z, x-homogeneous random potential V with 

(V(z,x)V(z',x.')) = 5(z-z') J e 4(x ~ x ' )p $(p)dp 
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where <3?(p) is the power spectral density and the product Vo$ in eq. stands for the Stratonovich 
product. Here and below the bracket (•) denotes the ensemble average w.r.t. the random medium. 
In view of the real-valuedness of V we may assume 

(2) $(-p) = $(p), Vp € R d 

Such a short-range-in- z potential arises in long-distance propagation when the longitudinal length 
L z of the wave beam is much larger than the correlation length of the random impurities in the 
medium resulting in rapidly fluctuating potential in the non-dimensionalized coordinates. And it 
is well known that in the short-range-correlation scaling limit a non-white-in-z multiplicative noise 
gives rise to a Stratonovich, instead of ltd, integral with respect to a white-in-z noise |12j . 

We are particularly interested in the problems of wave spread and singularity formation. To this 
end we shall use the phase-space approach of the Wigner distribution which is particularly useful 
in the regime of low Fresnel number 7 < 1 which can be viewed as a high-frequency limit. The 
main ingredient of our analysis is the phase-space variance identity by which we derive various 
estimates for the mean spread, defined in @, including precise behavior in the defocusing g < 
or critical case da = 2. To our knowledge, these are significant improvement over previous results 
(e.g. [Q, |14j . [3]) which are mostly for the linear problem. Wave spread is an important physical 
quantity in itself and in technological applications such as estimating the effect of nonlinearity on 
the information capacity of optical fiber communications [H], |1U| . 

Using the variance identity we further show that the scattering by the random potentials re- 
sults in conditions for singularity-with-positive-probability that are different from those for the 
homogeneous case. 

In the present setting the physical roles of t and z are interchanged: the variable t is space-like 
while the variable z is time-like and we will refer to it as "time" occasionally, especially in the 
discussion of finite-time singularity. 

2. Wigner distribution 
We consider the Wigner distribution of the form 

(3) W(z, x, p) = J e -«P** (z, x + ^) (z, x - ^) dy 

which is always real-valued and may be thought of as quasi-probability density on the phase space 
0,0. First its marginals are essentially the position and momentum densities 

Jw(z,x,p)dp = |*| 2 (z,x) = p(z,x) 

which are nonnegative; indeed any marginal on a ti-dimensional subspace is nonnegative everywhere. 
Also, the energy flux is given by 

(4) — (WT - = f pW(z,x,p)dp. 

2« jRd 

In the whole phase space, however, the Wigner distribution is not everywhere non-negative in 
general and hence can not be a true probability density. As 7 — > 0, the Wigner distribution has 
a nonnegative-measure- valued weak limit, called the Wigner measure [H]. Therefore the Wigner 
distribution is particularly useful for analyzing high frequency behaviors of waves. 
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3. The evolution equations 

From the NLS eq. Q it is straightforward to derive the closed-form (Wigner-Moyal) equation 
for the Wigner distribution in the ltd sense [I], [S] 

(5) d z W + p • X7 x Wdz + U^Wdz - Q^Wdz + V y Wdz = 0. 
Here the self-adjoint operator Q 7 is the Stratonovich correction term 

Q 7 ^(z,x,p) 

= J ^(q)7^ 2 [-2W(z, x, p) + W(z, x, p - 7 q) + W(z, x, p + 7 q)] dq 
and W 7 and V 7 are the nonlinear and linear Moyal operators, respectively 

llyW^Z, X, P) 

= i f e 1 ^- 1 [W(z, x, p + 7 q/2) - x, p - 7 q/2)] C/(z, q)dq, 
V 7 ^(z,x,p) 

= i y e 4q ' x 7 - 1 x, p + 7 q/2) - x, p - 7 q/2)] V(dz, dq) 

with 

[/ = <?p CT . 

Here we have somewhat abused the notation by writing U as the Fourier-transform-in-x of the 
function U(z,-) and V"(d-2r, dq) as the spectral measure of the white-noise potential V(z, •). The 
spectral measure t^(dz, dp) is related to the power spectral density <I> of V as 

'v(dz,dp)V*(dz',dei)\ = 5{z - z')S(p - q)$(p)dzdpdq. 

An important property of these integral operators is 

(6) J Q^Wdp = J U,,Wdp = j V^Wdp = 

which will be useful in deriving the energy law and the variance identity. 

A major advantage of formulating the wave process on the phase space using the Wigner distri- 
bution is that the (high-frequency) low-Fresnel number limit 7 can be easily obtained. Formally 
we see that as 7 — > 

U y W{z,x,p) ^U W{z,x,p) = V x C/(z,x) • V p VF(z,x,p) 

V 7 W(z,x,p) ^VoVF(z,x,p) = V x V(z,x) • V p lf( Z ,x,p) 

Q 7 W(z,x,p) ^Q W(z,x,p) =V p .DVpiy(z,x 1 p) 

with the diffusion coefficient 

(7) D = J $(p)p<g)pdp. 

We shall refer to eq. © as the nonlinear Wigner-Moyal -Ito (NWMI) equation when 7 > and as 
the nonlinear Liouville-Ito (NLI) equation when 7 = 0. 

Another advantage of working with eq. © is that one can use it to evolve the mixed-state initial 
condition. The mixed-state Wigner distribution is a convex combination of the pure-state Wigner 
distributions © described as follows. 
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Let {^ a } be a family of 1? functions parametrized by a which is weighted by a probability 
measure P(da) . Denote the pure-state Wigner distribution (j2J) by Wf^Q,]. A mixed-state Wigner 
distribution is given by 

(8) J W[* a ]P(da). 

The limiting set, as 7 — > 0, of the mixed state Wigner distributions constitute the nonnegative 
Wigner measures jS], [5]. The NWMI or NLI equation preserves the mixed-state form (JHJ). 



3.1. Basic properties. The NWMI and NLI equations conserve the total mass, i.e. 



dz 

and the 1? norm 



—N = 0, N = W(z, x, p)dxdp 



d_ 

dz 



J W 2 (z,x,p)dxdp = 0. 



One can absorb the effect of the total mass N into g by the obvious rescaling of W in the eq. [J5J). 
Henceforth we assume that N = 1. 

Let S x and S p be the spreads (or variances) of position and momentum, respectively 

(9) S x (z) = J \x\ 2 W(z,x,p)dxdp 

(10) S p (z) = J \p\ 2 W{z,x,p)dxdp. 

A natural space of initial data and solutions is the subspace W C L 2 (M. 2d ) with a finite Dirichlet 
form 

WQjWdxdp < 00 

and finite, positive (pre-ensemble-averaged) variances S x G (0, 00), S p € (0, 00). In addition, we 
shall also assume that the initial data have a finite Hamiltonian H £ (—00, 00) with 

(11) H = l -S p - -l-j J p° +1 dx, p(z, x) = J W(z, x, p)dp. 

The first term S p in H is the kinetic energy and the second term is the (self-interaction) potential 
energy. Note that in the presence of a random potential, the value of the Hamiltonian (|11|) is 
not conserved under the evolution (JSJ). Indeed, the ensemble-averaged value of the Hamiltonian 
increases in time, cf. ()29[) below. 
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4. Energy and variance identity 

Let us analyze the evolution of the mean Hamiltonian (H) and the average spread {S x ). First 
note the result of acting on the quadratic polynomials: 





(12) 


Zi-yX 


(13) 




(14) 


W 7 |x| 2 


(15) 




(16) 


w 7 lpl 2 



i J e iq x qC/(q)dq = V X U 







i J e iq ' x x • ql>(q)dq = x • V X C/ 
e iq x 2p • qt/(q)dq = 2p • V X U. 



A remarkable fact is that these results are independent of 7 > 0. 
Using (|12|) - (|16|) and © we obtain by integrating eq. © that 



■§^(S P ) = (^f V x U-2pWdxdp^ +j Q^\p\ 2 (W) dpdx 

d g 



. ha + i j ,>"■■ dx = {J YJ - V \Vdpdxy 
from which the evolution equation for the mean value of the Hamiltonian follows 

(17) A ( H ) = j Q 1 \p\ 2 (W)dpdx. 
We shall refer to ()17|) as the energy law. 

4.1. Variance identity. The variance identity has been long used to derive the wave collapse 
condition for the NLS in the homogeneous case [T^]. Below we reformulate it for the phase space 
evolution eq. ©. 

We have another set of remarkably simple properties of Q 7 independent of 7 > 0: 

(18) Q 7 x = 

(19) Q 7 p = y$(p-q)(p-q)dq = (by©) 

(20) Q 7 |x| 2 = 

(21) Q 7 x-p = x-y$(p-q)(p-q)dq = (by©) 

(22) Q 7 |p| 2 = - y ^(p-q) [|p| 2 - |q| 2 ] 

(23) = - J $(q)(2p-q)-qdq 

(24) = J $(q)|q| 2 (iq = 2i?. 

For the diffusion operator Qq we have R = 2 x tr(D) where the diffusion matrix D is given by ©. 
We shall use the above identities to perform integrating by parts in the derivation of the variance 
identity. 
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Combining the results from the previous section and Q18 |) -l|24 |) we obtain the rate of change of 

(S x ) 

(25) ^(S x ) = 2{S : 



where S xp is the cross-moment 



S X p = J x • pWdxdp. 
Differentiating S xp and taking expectation we obtain 

§~ Z (S* P ) = (S p )-f^J( P ^)dxd P 
Hence the second derivative of S x becomes 

(26) £ w . w+mzgnj^)^ 

Alternatively using the definition of H we can rewrite the variance identity (|2fi|) as 

(27) = Ma {H) + (2 - da) (S p ) 
which in the critical case da = 2 becomes 

(28) 1TT< 5 »> = 2 da(H). 

oz z 

Both (|26|) and (|27|) will be useful for estimating wave spread. 



5. Dispersion rate 

Although the medium is lossless, reflected in the fact that the total mass N = 1 is conserved, the 
value of the Hamiltonian, however, is not conserved by the evolution since the random scattering 
is not elastic due to time- varying nature of the random potential. Indeed, by Q17|) and (|24jl. the 
average Hamiltonian is an increasing function of time 

(29) £-{H) = R, (H)(z) = H(0)+Rz 

due to the diffusion-like spread in the momentum p. 

In the critical case da = 2, we obtain from eq. I)28|) and (|29|) the exact result 



(30) ^ (S x ) =4(H)= 4ff(0) + ARz 



dz 2 

before any singularity formation causes a possible breakdown of the variance identity. Integrating 
(|30j) twice we obtain the exact spread rate as stated below. 

Proposition 1. If da = 2 or g = 0, then 

(31) (S x ) (z) = S x (0) + 2S xp (0)z + 2H(0)z 2 + ^z 3 . 



The analogous result (S x ~ z 3 ) for the linear Schrodinger equation (d = 1, g = 0) with a random 
potential has been proved previously 

Next we consider the supercritical case and the defocusing case. We have from ()26(l that 

(32) £<4> = HH ) + ^kJ {p ^ )dx 



and hence 



d 2 

^ (Sx) < i(H) = 4iT(0) + 4i?z, for 5 (2 - da) < 
az 1 

d 2 

tt^ (Sx) > 4 m = 4HY0) + 4i?z, for 5 (2 - da) > 0. 



On the other hand, from ()27j) we obtain for any g 

d 2 

tt^ ^ x < 2da (H) , for 2 - da < 
d 2 

— (S*) > 2da (H) , for 2 - da > 0. 

Integrating the above inequalities twice, we obtain the following. 
Proposition 2. The following estimates hold 

(33) (S x (z)) < S x {0) + 2S xp (0)z + 2H(0)z 2 + ^Rz 3 , g(2 - da) < 

2 
3 

and 



(34) (Sr(*)> > S x (0) + 2S xp (0)z + 2H(0)z 2 + -Rz 3 , g(2 - da) > 



(35) (S^)) < S x (0) + 2S xp (0)z + daH(0)z 2 + ^i?z 3 , 2 < da 



(36) (S x (z)) > S x (0) + 25^(0)2 + daH(0)z 2 + — i?z 3 , 2 > da. 



da 
T 

Therefore 



Corollary 1. Assume g < (hence H > 0). Then 

(S x (z)) < S x (0) + 2S xp (0)z + (da V 2)H(0)z 2 + ^^Rz 3 

and 

(37) (S x (z)) > 5,(0) + 25^(0)2 + (da A 2)id(0)z 2 + ^-^Rz 3 . 

In other words the variance (S x ) is cubic-in-time in the defocusing case. In the subcritical case 
the cubic law is a lower bound while in the supercritical case the cubic law is an upper bound. 
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6. Finite-time singularity 



Finite-time singularity for the critical or supercritical (da > 2) self-focusing NLS equation in 
the homogeneous case is a well known effect [13 . In this case the singularity is the blow-up type 
S p , J |/o| <7+1 — ► oo. Here we call breakdown of the variance identity ()26 |) -()27 [) or the energy law 1)17(1 
as finite-time singularity and seek the sufficient conditions for singularity with positive probability. 
We show that in the supercritical case with additional assumptions the finite-time singularity is of 
the blow-up type in the sense that (S p ) , J (p CT+1 ) tend to infinity. As such the blow-up phenomenon 
discussed here is not necessarily a sure event but rather an event of a positive probability that S p 
and f p a+1 can exceed any fixed level in a finite time. 

For g > 0, da > 2 one can bound S x as in the inequality (|35|) 

(38) (S x (z)) < S x (0) + 2S xp (0)z + daH(0)z 2 + = F(z) 

and as motivated by the homogeneous case we look for the conditions when F(z) vanishes at a 
finite z > 0. 

A sufficient condition for F(z) to vanish at a finite positive z can be derived from that F(z) 
takes a non-positive value F(zq) < at its local minimum point zq > 0. The local minimum point 
zq is given by 



-H(0) + ^/H (0)2 - 2RS xp (0)/(da) 
(39) z = v — 



R 

Therefore we are led to the following conditions for singularity. 

Proposition 3. Assume da > 2,g > 0. The solutions of the NWMI or NLI equation develop 
singularities with positive probability at a finite time z* < zq given by under the condition 
F{ z o) < an d either one of the following conditions 

(40) S xp (0) < 

(41) S xp (0) > 0, H(0)<- ] f^5K. 

Remark 1. Clearly, the condition F(zq) < requires H(0) to be sufficiently below H by allowing 
the self-interaction energy 

--^ j p°+\0)dxdp 

to be sufficiently negative. 

Remark 2. Using, instead, inequality k3ty one can obtain an alternative expression 



(AO s -H(0) + y/g(0) 2 - RS xp (0) 
( 4 ^) z o = ^ 

and the corresponding conditions for singularity formation, namely either S xp (0) < or S xp (0) > 
0,H(0) < -y/RS xp (0). 

It may be the case that for a given initial condition the solutions develop singularity at different 
times depending on the realization of random potential. To analyze such an effect we need the 
probabilistic versions of variance identity and energy law which are much more involved. 
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6.1. Blow-up. We will follow the argument of |7j to show more explicitly the blow-up mechanism 
in the case with the supercritical, self-focusing nonlinearity and give a sharper bound on under 
certain circumstances. 

Proposition 4. Suppose da > 2,g > 0. Then under the conditions and 

2S (0) 

(43) (H(z*)) = H(0) + Rz, < 0, with z* = s ^ _ da) . 

d (S x ) jdz and (S p ) blow up before or at a finite time z*. 

Proof. Since blow-up is a local phenomenon, S x is a poor indicator of its occurrence. A more 
relevant object to consider is S p . From (|2T|) it follows that before z* 

(44) _ {Sx ) = 2— (S xp ) < (2 - da) (S p ) < 0. 

Since d (S x ) jdz = 2 (S xp ) < by (|H|) and 1)40(1 . d (S x ) jdz is a negative, decreasing function up to 
the time z#. Hence we have 

(45) < (S xp ) 2 < (S x ) (S p ) < S x (0) {S p ) , z E [0, z*} 
which implies 



(46) (S p ) > 



S x (0) 



Let A(z) = -(dS x (z)/dz) > 0,z < z*. We have from (jSJ), (@HJ) and ^ the differential 
inequality 

< 47 > i A> - CA '- c -%^)>° 



which yields the estimate 



A(0) 1 

A[z) - i-zCA(oy z< cmo)- z * 



and thus the blow-up of A(z) before or at z*. This along with (|46|) then implies the blow-up of S p 
at a finite time. □ 

The preceding argument demonstrates clearly the blow-up mechanism, namely the quadratic 
growth property (|4"T|) . 

Remark 3. It should be noted that, since H(zq) > 0, the condition implies that z* < zq, 

as given by either or j42\ ). Therefore, under the condition j43p , z* provides a sharper upper 
bound on the time of singularity than zq . 

Remark 4. Since (H(z)) is bounded over compact sets of z, Proposition^ implies the blow-up of 
the self-interaction energy f (p u+1 ^ in a finite time. 
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7. Conclusion 



We have presented an elementary approach for analyzing the nonlinear Schrodinger equation 
with a white-noise potential. We have focused on the ensemble-averaged quantities such as the 
variance identity and energy law in order to shed light on two problems: the rate of spread and the 
singularity formation. 

We have shown that the ensemble-averaged spread in the critical or defocusing case follows the 
cubic-in-time law while in the supercritical and subcritical focusing cases the cubic law becomes 
an upper and lower bounds respectively. In a separate publication we will use these estimates 
to analyze the limitation on channel capacity in optical fibers due to self-phase and cross-phase 
modulations [S]. 

We have also found singularity conditions in the critical and supercritical focusing cases. And we 
show the finite-time singularity in the supercritical case is of the blow-up type. The singularity /blow- 
up discussed in the present paper is not necessarily a sure event but that of a positive probability. 

References 

[1] L. Bunimovich, H.R. Jauslin, J.L. Lebowitz, A. Pellegrinotti and P. Nielaba, Diffusive energy growth in classical 

and quantum driven oscillators. J. Stat. Phys. 62(1991), 793-817. 
[2] A. Debussche and L. Di Menza, Numerical simulation of focusing stochastic nonlinear Schrodinger equations. 

Phys. D 162 (2002), 131-154. 

[3] E.B. Erdogan, R. Killip and W. Schlag, Energy growth in Schrddinger's equation with Markovian driving. 

Commun. Math. P/ij/s.240(2003), 1-29. 
[4] A. Fannjiang, White-noise and geometrical optics limits of Wigner-Moyal equation for wave beams in turbulent 

media. Comm. Math. Phys.254:2 (2005), 289-322. 
[5] A. Fannjiang, Self-averaging in scaling limits for random parabolic waves. Arch. Rat. Mech. Anal. 175:3 (2005), 

343 - 387. 

[6] P. Gerard, P.A. Markowich, N.J. Mauser and F. Poupaud, Homogenization limits and Wigner transforms, Comm. 

Pure Appl. Math. I (1997), 323-379. 
[7] R.T. Glassey, On the blow-up of solutions to the Cauchy problem for nonlinear Schrodinger equations, J. Math. 

Phys. 18(1977), 1794-1797. 
[8] L. Kazovsky, S. Benedetto, A. Willner, Optical Fiber Communication Systems, Artech House, 1996. 
[9] P.L. Lions and T. Paul, Sur Jes mesures de Wigner. Rev. Mat. Iberoamericana 9(1993), 553-618. 
[10] P.P. Mitra and J.B. Stark, Nonlinear limits to the information capacity of optical fibre communications, Nature 

411 (2001), 1027-1030. 

[11] A. C. Newell and J. V. Moloney, Nonlinear Optics, Addison-Wesley, Redwood City, CA., 1992. 
[12] G. Papanicolaou and W. Kohler, Comm. Pure Appl. Math. 27 (1974), 641-668. 

[13] C. Sulem and P.-L. Sulem, The Nonlinear Schrodinger Equation: Self-Focusing and Wave Collapse. Springer, 
New York, 1999. 

[14] S. Tcheremchantsev, Transport properties of Markovian Anderson model, Commun. Math. Phys.l9&:l (1998), 
105-131. 

E-mail address: cafannjiang@ucdavis.edu 

Department of Mathematics, University of California, Davis 95616-8633 



10 



